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UITTREKSEL 



Aan de hand van een reeks xeflectiemetingen aan een vloeistof worden 
parameters berekend die een deeUgesgrootteverdeling van deeltges in de 
vloeistof beschrijven. In elke reflectiemeting wordt een signaalbundel in de 
vloeistof wordt opgewekt en een waarde van een eigensdiap zoals de 
amplitude van een reflectie aan een deelije in de signaalbimdel gemeten. 
Met behulp van een maximum likelihood schatting worden de parameters 
van de deellgesgrootteverdeling geschat, op basis van een uitdrukking voor 
een waarschijnlijkheid van de gemeten waardes als fimctie van de gemeten 
wa^irdes. De gebnufcte uitdrukking bevat een eerste factor voor de 
waarschjjnlijkheid van een reflectiemeting waarvan een reflectie met de 
gemeten waarde onderdeel van uitmaakt, gecorri^erd met een tweede 
factor voor de waarschijnlijkheid dat niet tevens een reflectie met een 
dominerende waarde van de eU^nsdhap, die de gemeten waarde zou 
maskeren, deel iiitmaakt van de reflectiemeting. 
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Titel: Werkwijze en ixirichfcixig voor het verkrijgen van informatie over de 
grootteverdeling van macroscopische deeltjes in een vloeistof 

De uitvinding betreft een werkwijze en inrichting waarin, met bebulp van 
weerkaatsing van een vdtrasone pulsen door macroscopische deeltges zoals 
van oliedruppels of zandkorrels in water, informatie over de 
grootteverdeling van de deellges in een vloeistof bepaald wordt. 

5 

Europese octrooiaanvrage No. 0801305 beschrijffc een inrichting 
voor het karakteriseren van een suspensie. De inrichting wekt xdtrasone 
pidsen op in eenbundel in een vbeistol Reflectie van een puis aan een 
individueel deellje in de vloeistof resulteert in een echo. Als het deellge zich 

10 op een gegeven positie bevindt is de ampUttwie van de reflectie e6n op 64n 
aan de grootte van het deellrje gekoppeld. Bij metingen aan reflecties van 
achtereenvolgende pxdsen op deeltjes van verschillende grootte worden 
zodoende verschillende amplitudes gemeten. 

De bekende inrichting maakt een histogram van tellingen van het 

15 aantal keren dat reflecties met verachillende amplitude voorkomen. Het 
histogram bevat informatie over de deeltjesgrootteverdeling, dat wil zeggen 
de concentratie van deeUjes als functie van de grootte van de deelljes. Met 
informatie over de deelljesgrootteverdeling op zqn beurt kan bgvoorbeeld de 
totale concentratie van deelfges in de vloeistof bepaald worden. 

20 Europese octrooiaanvrage No. 0801305 beschrjjffc een techniek cm 

informatie over de grootte uit de verdeling van de amplitudes te halen. Er 
wordt een tdtdrukking gebruikt voor de waarschijnlijkheid dat reflecties 
met verschillende amplitudes optreden. Deze uitdrukking relateert de 
waarschijnlijkheid aan de deeltjesgrootteverdeling. Als de reflecterende 

25 deeltjes zich aMjd op dezelfde plaats in de bundel zouden bevinden dan zou 
het aantal reflecties met een bepaalde amplitude evenredig zijn met de 
fractie aan deeltjes met de grootte die tot die amplitude leidt. 



024984 



V 



De informatde over de deeltjesgrootte verdeling zit echter 
versmeerd over de amplitudes, omdat de refiectieamplitude van een deeUje 
behalve van de grootte van het deeltjes ook afhankelj|k is van de positie van 
het deellje in de bundeL Hoe verder het deeltje zieh bij reflectie nit het 
5 centrum van de bimdel bevindt, hoe kleiner de amplitude. De gebruikte 
uitdrukking voor de waarschijnlgkheid van een reflectiemeting met een 
bepaalde amplitude versmeert de deeltjesgrootteverdeling over de 
amplitudeverdeling om dit efiect tot uitdrukking te brengen. 

Met behulp van een maximum likelihood techniek worden 

10 parameters van de deeltjesgrootte verdeling geschat die de zo uitgedrukte 
waarschqnlilkheid van de daadwerkelQk gemeten aantallen deellges 
maximaUseren. Zodoende wordt de versmerxng ongedaan gemaakt. In 
concreto worden de parameters zo gekozen dat een som geminimaUseerd 
wordt van de kwadraten van de verschillen tussen gemeten aantallen 

16 reflecties met amplitudes in verschiUende amplitudebereiken en voorspelde 
aantallen. Een dergel^jke som is een indicatie voor de waarschonlilUieid 
van de combinatie van gemeten reflecties, maar uiteraard kunnen ook 
andere indicaties van de waarschijnUjkheid gebruikt worden, zoals het 
product van de waarschijnlijkheden van de verschiUende gemeten 

20 amplitudes. 

E!en noodzakelijk voorwaarde voor de gebruikte techniek is dat de 
reflecties van individuele deellges onderscheiden kunnen worden. Als 
reflecties. van meerdere deelljes door elkaar gemeten worden, dan raakt een 
deel van de reflecties gemaskeerd. Zo kan reflectie van een groter deeltje 

25 reflectie van een kleiner deelfje onzichtbaar maken, maar ook kan reflectie 
van een kleiner deeltge reflectie van een groter deeltje maskeren als het 
grotere deeltje zich voldoende veel verder vdt het centrum van de bundel 
bevindt dan het kleinere deehrje. 

Zodoende geeft de techniek alleen brj voldoende lage concentraties 

30 befxouwbare resultaten. De maximaal bruikbare concentratie kan verhoogd 
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worden door hot minimaliseren van hot voltune waarin deeltjes 
ononderscheidbaar gemeten worden, bqvoorbeeld door hot gebruik van een 
gefocusseerde bundel en Meine tijdswindows waarin edio*8 geaccepteerd 
warden. Er zijn echter grenzen aan de toepasbaarheid van dergelijke 
5 technieken. Daarom beperkt de eis van gescheiden waameming van 
individuele deeltjes de toepasbaarheid van de techniek. 

Het is een doel van de nitvinding om het gebruik van tiltrasone 
metingen voor de karakterisatie van deeWrjesgrootteverdeling in vloeistoflfen 

10 bruikbaar te maken voor bogere concentraties. 

De huidige tiitvinding berust op de realisatie. dat bet mogelgk is om 
met een statistisch model bet efifect van bet missen van deelljes te 
corrigeren en dat daardoor de bekende meettechniek ook voor grotere 
concentraties gebruikt kan worden. 

15 De uitvindiz^; voozziet in een werkwqze volgens de condusie 1. In 

de werkwijze worden de parameters van de deeltjesgrootte verdeUng 
geschat met een maximum likelihood techniek, waarbij gebruik gemaakt 
wordt van een uitdrukking voor de waarscbijnlijkheid van meting van 
reflectieamplitudes in termen van twee fectoren. De eerste fiactor drukt de 

20 waarscbijnlijkheid Po(A) uit van een reflecliemeting waarvan een reflectie 
met de waairde A onderdeel van uitmaakt, onafhankelijk van de vraag of 
deze waarde gemaskeerd wordt doordat tevens een reflectie met een grotere 
waarde opteeedt. Deze eerate fiactor kom in wezen overeen met de hole in de 
stand der tecfaniek gebruikte uitdrukking voor de waarscbijnlijkheid. 

25 De tweede factor in de volgens de uitvinding gebruikte uitdrukking 

omvat de waarscbijnlijkheid dat niet tevens een reflectie met een 
dominerende waarde van de eigenschap, die de gemeten waarde zou 
maskeren, deel uitmaakt van de reflectiemeting. Zodoende wordt het effect 
van maskering hij de schatting van de parameters meegenomen. Bij de 
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likelihood schattiiig woxdt zodoende zowel het effect van 
versmerix^ als het effect van maskering samen meegenomen. 

Deze en andere doelen en voordelige aapecten van de uitvinding zullen 
5 worden beschxeven aan de hand van de volgende figuren 

Figuur Itoont een inxichting voor het karakteriseren van een vloeistof 
Figuur 2 toont aantallen xeflecties met verschillende amplitudes 
Figuur Stoont een gereconstrueerde deeltjesgrootteverdeling 
10 Figuxir 4toont een berekende concentratie 

Figuur Stoont een flow chart van een werkwijze om parameters van een 
deeltjesgrootteverdeling te schatten 

F^ur 1 toont een inxichting voor het karakteriseren van een vloeistof, De 
15 inriclitdng omvat een vloeisto£kanaal 11, conteole eenheid 10, een 

si^naalgenerator 12, een ultrasone transducer 14, een detector 16 en een 
rekeneenheid 18. Ck>ntrole eenheid 10 is aan signaalgenerator 12, detector 
16 en rekeneenheid 18 gekoppeld. Signaalgenerator 12 heeft een uitgang die 
aan transducer 14 gekoppeld is. Detector 16 heeft een ingang die aan 
20 transducer 14 gekoppeld is. Detector 16 heeft een uitgang die aan 

rekeneenheid 18 gekoppeld is. Transducer 14 is aan kanaal 11 gekoppeld 
voor het opwekken van een ultrasone bimdeL 

In bedxijf stroomt een vloeistof met daarin zwevende deeltjes door kanaal 
11. Ciontrole eenheid 10 tr^gert achtereenvolgende reflectiemetingen in elk 

26 waarvan signaalgenerator 12 een puis genereert die door transducer 14 in 
een ultrasone golf in kanaal 11 wordt omgezet. Transducer 14 vangt 
reflecties van de puis pp en voert resulterende si^alen aan detector 16 toe. 
Detector 11 selecteert signalen die in een bepaald tijdswindow ten opzichte 
van de opgewekt puis optreden. Als gevolg van de geluidssnelheid in 

30 vloeistof definieert het tijdswindow een zone vanwaaruit deeltjes aanleiding 
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Inmnen geven tot reflecties. Detector 16 meet de amplitude van het 
xeflectiesigiiaal (indien aanwezig) in het tgdswindow en voert de gemeten 
waarde van de amplitude toe aan lekeneenheid 18. 
Bekeneenheid 18 verzamelt histogram informatie van de amplitudes* 
5 Hiertoe gebruikt rekeneeuheid een aanlal amplituderas^s en accumuleert 
rekeneenheid tellingen van het aantal reflectiemetingen die 
reflectieamplitudes in de respectieve raises hebben opgeleverd. Nadat 
zodoende een groot aantal reflectiemetingen is uitgevoerd berekent 
rekeneenheid 18 uit de histogram informatie verdere informatie over de 

10 deeltjesgrootte verdelii^. 

F^guur 2 toont een gesimuleerd voorbeeld van histogrammen voor 
vloeistofiEen waarin verschillende ooncentraties van deelljes aanwezig zqn, 
terwgl die deeUjesgrootteverdeling foCD) telkens hetzelfde is (dat wil zeggen 
dat de fractie van de deelljes met een grootte tussen D en D+dD telkens 

15 fi)CE>)dD is voor elke waarde van D). Zichtbaar is dat» behalve een toename 
van de aantallen reflecties ook de vorm van de amplitude verdeling 
verandert biij toenemende concentratie. 

Figuur 3 toont met gestippelde lijnen schattingen van de 
deelljesgrootte verdeling die uit de amplitude histogrammen kan worden 

20 berekend op basis van de bekende techniek, voor een aantal verschillende 
concentraties. Tevens wordt de werkelijke verdeUng getoond. De bekende 
tedbniek houdt geen rekening met het effect van de concentratie op de vorm 
van de amphtudeverdeling. Het residtaat is dat de geschatte 
deeltjesgrootteverdelingen onderling ve»chillen en bij hogere ooncentraties 

25 afwijken van de echte deelfjei^rootteverdeling. 

Da uitvinding voorziet in een verbetering van de bekende techniek 
om de deeltjesgrootteverdeling te schatten, zodat ook voor het ejBEect van de 
concentratie van de deeltjes gecorrigeerd kan worden. Deze verbetering 
berust op de realisatie dat de afwijking van de vorm van de histogrammen 

30 van figuur 2 een gevolg is van maskering van reflectieamplitudes die bij de 
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een zeflectiemBtixig wel optreden, maar in de telling gemist wordon omdat 
Wj dezelfde reflectiemeting 6en of meer reflecties met een giotere amplitude 
optzeden. 

Daartcw wordt een nitdrukkii^ gebroik vror de 
5 waaxschijii]q]Eheidsverdeling P(A.)dA dat een reflectiemeting een gemeten 
amplitude oplevert die in een infinitesimaal interval tuBsen A en A+dA ligt. 
Deze idtdrukking heeft de vorm van 

P(^=Po(A)Q(A) 

10 

Hierin ooriespondeert Po(^ met de ongemaskeerde 
waar8db:rjnlqkheidsverdeling, dat wil zeggen waarschijnljjkheid van 
aanwesdgheid van een leflectie, los van de vraag of deze gemaskeerd wordt 
door een andeie reflectie met een grotere amplitude. De factor Q(A) drukt de 

15 waarschijnigkheid uit dat er in de dezelMe reflectiemeting geen deelljes 
reflecties met grotere of gelijke amplitude dan A verooraaken. 

De &ctor Po(A) komt overeen met de waarschijnlijkheidsverdeling 
die in de stand der techniek gebruikt wordt voor de gebele 
waarsch^jnlijkheid verdeling P(A) dat een amplitude A in het interval A tot 

20 A+dA ook daadwerkelijk gemeten wordt (dwz onder verwaarlozing van de 
mogelijkheid dat een reflectie gemaskeerd wordt). Deze factor Po(A) kan op 
op zichbefcende wqzse in de deellgesgnmtteverdeling fi>(D) worden 
tiitgedrukt. Een tiitdrukkii^ voor deze &ctor is bijvoorbeeld 

25 Po CA)= C JdD foOD) G(A| D) VmeasCD) 

Ifierin is C de concentratie van de de^es (het gemiddeld aantal deelljes 
per volume-eenheid). G(A| D) is de oonditionele waarsdiijnlijkbeid dat er 
een reflectie met amplitude A detecteexbaar is, als een deeltje met een 
30 grootte D een detecteerbare reflectie oplevert. De conditionele 
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waarschqnlijkhfiid G(A| D) drukt hot efEect uit dat een deelije met een 
grootte D op versGhiUende afstanden van het centrum van de bundel 
aanleiding geeft tot verscbfllende reflectie amplitudes. Vm^W) is het totale 
vloeistofvolume vanwaaruit een deelige met grootte D een detectee Aare 
5 reflectie kan opleveren. Dit volume wordt in de diepte gedefinieerd door het 
tgdswindow vanwaaruit tijdens de meting reflecties geacoepteerd worden en 
in de breedte door de afval van de bundelintensiteit tot een punt waaxop de 
amplitude van de opgewekte reflecties niet meer boven het ruisniveau An 
uitsteekt. 

10 In het veel voorkomende geval van een Gaussisch verdeelde bundel 

kan het volume VmeasCD) worden uitgedrukt als 



Vnieas(D)=Il AZ l0g(Ao(D)/ AayZK 

15 Hierin is Az het diepteinterval waaruit reflecties geaocepteerd worden, k is 
de standaard deviatie die bepaald wordt door de snelheid waarmee de 
intensiteit van de bundel als functie van de afetand tot het centrum van de 
bundel afvalt. Ao(D) is de maximale amplitude die een reflectie van een 
deellge met grootte D zou opleveren als het deeltje zich in het centrum van 

20 de bundel zou bevinden. 

Voor een Gaussische bundel kan ook een uitdrukking voor de 

conditionele waarsch^nlijkheid G(A. I D) worden gegeven: 



G(A|D)= l/{A*log(Ao(D)/AB)> 

Deze uitdrukking geldt mits Ao(D)>A>An . Voor amplitudes groter dan de 
maximaal met een deeltjesgrootte D haalbare amplitude Ao(D) of Heiner dan 
het ruisniveau An is de conditionele kans op detectie G(A | D) gelijk nuL 
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Het zal ovexigem duidelgk zijn dat vocw btindels met een andexe 
dan Gaussische vexdeling ovexeenkomsti^ andere uitdrukMngen voor 
VmeasCD) en GKA I D) gebruilrt zuUen woiden. 

De meuw gelntxoduceerde fsctox Q(A), die het effect van maskezmg 
5 van lefiecties door andexe xeflecdes met een gxotere of gelijke amplitude dan 
de amplitude A uitdrukt, kan eveneens woxden uitgedrukt in texmen van de 
concentratie C en de deelijesgrootteverdeling &(D). Maatgevend hiervoor is 
de integraaluitdriikking 

10 W = JdA'/dD &(D) G(A' I D) 

Hierin woxdt de integraal over de amplitude genomen van het ruisniveau tot 
aan de amplitude A, waaxvoox de kans op detedae beiekend wordt. De &ctor 
Q kan worden benaderd door 

15 

Q(A)=exp[-C(l.W(A))] 
Of bij benadexir^ door een £actor 

Q C/^ = [ 1-p +p (l-W(A)Armeas) ]» 

20 

Waaxin de som over k genomen wordt van 1 tot n; n is een geheel getal dat 
niet gfot&t is dan het mfi-Hmaal aantal deellges dat Igrsisch past in het 
volume vanwaaruit meetbare leflecties optreden, maar ruimschoots groter 
dan het gemiddeld aantal deelljes dat meetbare reflecties oplevert en 
25 p=CVmeas/n, waarbij Vmeas het volume is van waaruit de grootste deeltges 
nog echos boven het ruisniveau kunnen opwekken. - 
Een altematieve xiitdrukkii^ voor de P o (A) factor is bajvoorbeeld 

Po (A)= C JdD fo(D) cfV(AD)/ dA 

30 
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Hierin is C de ooncentratie van de deelljea (het gemiddeld aantal deelijes 
per volume-eenheid). V(A,D) ia het volume waaruit een deeltje met grootte D 
een echo met een amplitude groter of gelijk A Qplevert. Dit vdlmne wordt in 
de diepte gedefinieerd door het tqdswmdow vanwaaruit tqdens de meting 
5 refiecties geaccepteerd worden en in de breedte door de afiral van de 
bundelintensiteit met de afetand tot het centrum van de bundeL 

In het veel voorkomende geval van een Gaussisch verdeelde bundel 
kan het volume V(A,D) worden uitgedrukt als 

10 V(A, D)=it Az log(Ao(D)/ A)/2k 

Hierin is Az het diepteinterval waaruit refiecties geaccepteerd worden, k is 
een parameter die de snelheid karakfceriseert waarmee de intensiteit van de 
bundel als functie van de a&tand tot het centrum van de bundel afvalt. 

15 Ao(D) is de maxunale amplitude die een reflectie van een deeltie met grootte 
D zou opleveren als het deeltje zich in het centrum van de bundel zou 
bevinden. Het zal overigens duidelijk zijn dat voor bundels met een andere 
dan Gaussische verdeling overeenkomstig andere uitdrukkingen voor V(A, 
D) gebruikt zal worden. 

20 De nieuw geinteoduceerde fector Q(A), drukt de waarschijnl^jkheid 

uit dat er naast de echo met amplitude A geen echo's met grotere ampUtude 
optreden. Dit beschrqft het efiBsct van maskering van refiecties door andere 
refiecties met een grotere of gelijke amplitude dan de amplitude A. De &ctor 
Q(A), If"-" eveneens worden uil^edrukt in termen van de ooncentratie C en 

25 de deehjesgrootteverdeling fbCD). Maatgevend hiervoor is de 
integraaluitdrukking 

WCA)= /dD&CD)V(A.D) 
30 De £Bu:tor Q kan worden benaderd door 
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QCA)=ejq»[-CW<A)] 

Aan de hand van de tdtdxukkmg P(A.)= Po(^ Q voor de 
5 waarachgiiiiildieidsverdelixig van metii^n van refledaes met amplitude A 
die niet gemaskeerd woiden door reflecties met groteze amplitude woidt uit 
het stelsel van lefiectie amplitudes dat in een gioot aantal nratingen 
gemeten is de deelljes grootte verdeling fo(D). Dit gebeurt bijvoorbeeld door 
een schatting fo(D) te 2»eken die de som E minimaliseert van de kwadraten 
10 tussen de voorspelde aantallen metingen Npted van amplitudes in een reeks 
amplitude intervallen lAi en de daadwerkelijk gemeten aantallen Nmeas 

E - B (NnmasOAO- NpiedOAO)^ 

15 Deze som is in feite een indicaMe voor de waaxschijnlqkheid van het 

gemeten stelsel van amplitudes, zodat de deeligea grootte verdeling fo(D) 
gekozen wordt die dit stelsel het meest waarschijnlijk nmakt. Uiteraard kan 
in plaats van de som van de kwadraten ook met behtdp van andere 
indicalieB van de waarschqnlijkheid gemaximaliseerd worden, bijvoorbeeld 

20 met behulp van een product 

niPo(Ai)Q(A«) 

van de "voaarschgnlijkheden van de ve»chillende metingen At. 
25 Bg deze schattii^ wordt bi^rvoorkeiir een parametrisatie van de 

deeltgesgrootteverdeling fi>CD) gebruikt, bijvoorbeeld in de vorm van een 
multimodal verdeling 

fi>(D) = 2i Ci e3q)(-(D-Di)2/2<yi2 

30 
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(i=l, In dit geval worden de parametexs Ci Di en zo geschat dat de 
volsens de gel»ruikte indicatie de waarschqnlqkheid van de hst gemeten 
stelsel refledae amplitudes gemaxunaliBeerd wordt. Technieken voor 
dergelgke maximalisaties zijn in hun algemeenheid bekend, en door deze 
5 toe te passen op een nitdmkkmg Pa(A) Q (A) voor de waarschijnUjkheid van 
gemeten (en niet gemaskeerde) reflectieamplitades verkrijgt men een 
schattingvoor de parameters. 

Figuur 3 toont tevens verschillende deellgesgrootteverdelii^n die 
zodoende nit de aantallen reflecties in de voorbeelden van figuur 2 berekend 
10 worden. verdelingen vallen binnen de nauwkeurigheid van de figutur 
nagenoeg samen met de werkeUjke verdeling. Het zal duidelijk zijn dat 
2X)doende een betere schatting van de deelijesgrootteverdeling verkregen 
wracdt. 

Figuur 4 toont de gevolgen van de schatting van de deelljesgrootte 

15 verdeling voor de conoentratie. De ooncentratie is de firactie van de vk)eistof 
die door de deeltjes wordt gevormd. Horizontaal staat de werkelijke 
concentratie en verticaal de uit de (gesimuleerde) metingen berekende 
ooncentratie. De getrokken lijn geeft aan waar berekende en werkelijke 
concentratie gelijk zijn. De drkels geven de met de nieuwe techniek 

20 berekende conoentratie en de krui^es de met de stand der techniek 

berekende concentraties. Uit de figuur blijkt dat met name voor hogere 
ooncentraties duidel:Qk betere resultaten verkregen worden. 

Het zal duidelqk zijn dat de uitvinding niet gebonden is aan deze 
spedfieke wijze van schatten van de deeltjesgrootteverdelii^. Zo kan 

25 bijvoorbeeld een andere vorm van parametriseren van de deeUjesgrootte 
verdeli]]^ gebruikt worden. bjjvoorbedd aan de hand van de waarden (en/of 
afgeleiden) van de deeUjesgroottev^rdeling voor een aantal deeUjesgroottes 
en een interpolatie van de deeltjesgrootteverdeUng tussen deze 
deel^esgroottes, of bijvoorbeeld in de vorm van een quotient van twee 

30 polynoxnen, waarvan de coSfiEidenten de parameters vormen. Ook kunnen 
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benaderingen voor d© versdiillende termen in de wiskundige formule voor 
de kans P(LA^ gebndkt woxden. 

Vetder kan bgyoorbeeld eerst een schatting gemaakt worden van 
de fector Q(A) als functie van de ampUtude A (bijvoorbeeld door provisorisch 
5 de deelijesgrootte verdeling te schatten en vervolgens de fector met die 
provisorisch geschatte verdeling nit te rekenen) en vervolgens kan men de 
gemeten aantallen reflecties NmeasCA) met verschfllende amplitudes 
corrigeren met deze fector Ncorr(A)= NmeasCAVQCA). Uit de gecorrigeerde 
aantallen kan vervolgens op de vat de stand der tedmiek bekende wijze de 
10 deelijesgiootte verdelxi^ geschat worden. 

Hoewel de nitvinding beschreven is in termen van de e^lidete 
deeWjesgrootteverdeling fi)(D) zal bet dnidelijk zijn dat bet voor bepaaMe 
toepassingen niet nodig is om deze deeHgesgrootte verdeling ook e^lidet nil 
te rekenen. Wanneer alleen een grootheid zoals bijvoorbeeld de 
15 massadichtbeid aan deelljes nodig is, die volgt uit de integraal 



idD fo(D) m(D) 

dan kan deze grootheid, die een parameter van de deelldesgrootte verdeling 
20 is ook implidet berekend worden bij bet schatten van de parameters van de 

deelljesgrootteverdeling die het stelsel van gemeten reflectiB amplitudes het 

meest waarschijnlqk maakt. 

De gebruikte techniek is niet beperkt tot metingen aan amplitudes, 

maar kan uitgebreid worden op elke vorm van metingen van eigenschappen 
25 van de reflecties waarb\i, btj aanwezii^ieid van meerdere deelljes die 

reflecties veroorzaken, het effect van een dominerend deellje het effect van 

de overige deelljes maskeert. 

De techniek is met name ook toepasbaar op vloeistoffen die 

verschillende soorten deeltjes bevatten, wanneer onderscheid gemaakt kan 
30 worden van welk sooit deellje tot een reflectiemeting aanleiding geefl;. Bit 
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ondeiseheid kan bqvoorbeeld aan het zeflectieaignaal zeK zichibaax zijn, 
bjjvooxbeeld door het token van de eexste piek in een ontvangen ultrasoon 
reflecfciesignaal in le^onse op een ultrasone puis met een bepaald teken. In 
het geval van water waaxin oliednqn>el8 en zandkorrels zijn (vgelost, 
5 bijvoorbeeld, kunnen lefiecties aan zandkonels en oliednflppds 
onderscheiden worden aan het teken van de reflectie. 

De gecorrigeerde waarsdbijnlijkheidsverdeling Pj(A)dA voor een 
meting van een reflectie van een deellje van type j, met een amplitude in 
een bereik tussen A en A+dA kan in dit geval beschreven worden met 

10 

Pj^Pjp(A)niQiCA) 

Hierin beschrijven verschillende factoren Qi (A) voor elk verschillend soort 
deellaes i de waarschijnlijkheid dat geen deeltjes van het soort i tot detectie 

15 leiden met grotere amplitude dan A. De factoren Q i (A) worden elk op de 
bovenbeschreven manier uitgedrukt in de deeHjesgrootteverdeling fi)i(D) va 
een respectievelqke soort deeltjes i, evenals de waarschijnlijkheden Pjo(A). 
Met deze uitdrukkingen voor Pj(A) voor verschillende soorten deeltges en ee 
stelsel metingen van xeflectieamplitudes die geklassificeerd zijn naar het 

20 betrokken soort deeUjes kunnen op de bovenbeschreven wijze parameters 
van de deeHjesgrootteverdehngen fi)i(P) geschat worden. Zo kunnen 
bijvoorbeeld deellje^rootteverdelingen geschat worden die een 
samengestelde som minimahseren van kwadxaten van aftirijkii^n 

25 E = Si In. (NimeasOAn)- Ni,«ed(IABi))2 

van de aantallen Nimeas van gemeten reflecties van deeltjes van soort m in 
verschillende amplitudebereiken lAm en de voorspelde aantallen Nipred 
gemeten reflecties in deze bereiken. 
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Biguur 5 toont samenvattend de werkwijze voor het echatten van 
de deelljesgrootteverdeling. In een eerst stap 61 wordt een ptds opgewekfc en 
in een tweede sitap 52 wordt gedetecteerd of er in een bepaald tijdswindow 
na de puis een leflecfcie boven het niisniveau terugkomt. In een derde stap 
5 53 wordt bepaald in welk van een aantal iatervallen lA de amplitude valt en 
wordt in een geheugen een telling voor het betrokken interval opgehoogd. 
Vanuit een vierde stap 54 worden de eerdere stappen herhaald totdat een 
bepaald aantal N reflectie espeiimenten is uitgevoerd. Vervolgens wordt een 
vijfde stap 56 uitgevoerd waarin parameters van een initiele 

10 deelio^sgrootteverdeling gekozen wordt. Hiermee wordt in een zesde stap 56 
met behulp van de formula P(A)=Po(A) Q (A) een voorspelling van de 
gemeten aantallen uitgerekend en wordt een verschil met de gemeten 
aantallen bepaald. Als het verschil voldoende klein is stopt de werkwgze in 
een zevende stap 57. Anders worden in een achtste stap 68 de parameters 
. 15 aangepast, zo dat het verschil naar verwachting zal verminderen en wordt 
de werkwijze vanaf de zesde stap 56 herhaald. Uiteraard wordt schatting in 
de werkwijze bqvoorkeur met een computer uitgevoerd. 
Het zal duideHjk zijn dat de uitvinding niet is beperkt tot de beschreven 
werkwijze, zo hoeft bijvoorbeeld in de zesde stap 56 de het verschil niet 

20 expHdet uitgerekend te worden, maar kan in plaats daarvan infi^naatie 
uitgerekend worden die nodig is om de aanpassing van de parameters te 
kiezen en kan de werkwqze stoppen in de zevende stap 57 als de aanpassii^ 
Meiner dan een drempel is. 
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CONCLUSIES 

1. Werkwqze voor het meten van parameters die een 
deeltjesgrootteverdeling van deelljes in een vloeistof beschrijven, welke 
werkwqze de stappen omvat van 

- het uitvoeren van een reeks reflectiemetingen, in elk waarvan een 

6 signaalbimdel in de vloeistof wordt opgewekt en een waarde (A) van een 
eigenschap van een reflectie aan een deeltje in de signaalbimdel wordt 
gemeten; 

- het uitvoeren van een maximum likelihood schatting van de parameters 
gezien een combinatie van de gemeten waardes (A), op basis van een 

10 uitdrukking voor een waarsduonhjldieid van de gemeten waardes als functie 
van de gemeten waardes, welke uitdrukking een eerste fiactor (P) bevat voor 
de waarschijnlijkheid van een reflectiemeting waarvan een reflectie met de 
gemeten waarde onderdeel van uitmaakt» gecorrigeerd met een tweede 
fsLCtor (Q) voor de waarschgnljjkheid dat met tevens een reflectie met een 

16 dominerende waarde van de eigenschap, die de gemeten waarde zou 
maskeren, deel uitmaakt van de reflectiemeting. 

2. Werkwijze volgens conclusie 1, waarin de eerste factor (P) de 
deelljesgrootte verdeling omvat, versmeerd met een kansverdeling dat een 
deelfge met een bepaalde grootte leidt tot een reflectiemeting waarvan een 

20 reflectie met de gemeten waarde onderdeel van uiianaakt. 

3. Werkwijze volgens conclusie 1 of 2, waarin de tweede factor (Q) de 
waarsdugnlgkheid dat een reflectie met een andere dan de gemeten waarde 
onderdeel van idtmaakt van een reflectiemeting, geintegreerd over een 
bereik aan andere dan de gemeten waardes. 

25 4. Werkwijze volgens conclusie 3, waarin de tweede factor nagenoeg 

overeenkomt met exp(- C JdA^JdD &CD) G(A* I D)), waarin D een 
deeltjesgrootte is, C een concentratie van de deelljes, fo(D) een dichtheid 
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van deellges met deeltjesgrootte D en G<A | D) is de conditionele 
waarschijnlijkheid dat er een reflectie met amplitude A detecteerbaar is» als 
een deellge met een grootte D een detecteerbare reflectie oplevert. 

5. Werkwijze volg^ns 6en der voora%aande oondusies, waaxin de 
5 ywa-riTtiiiTn likelihood schattii^; het maken van tellingen omvat van 

aantallen reflectaemetingen waaxin de gemeten waarde in lespectieve 
waardeintervallen vallen, en de schatting zo gekozen wordt dat een 
samenstel van afwijkingen tussen de tellixi^en in de verschiUende 
intervallen volgens de waarschijxilijklieid als fimctie van de gemeten 
10 waardes voorspelde tellingen geminimaliseerd wordt. 

6. Werkwijze volgens een der voorafgaande condnsies, waarin bij de 
reflectiemetingen onderscheid gemaakt wordt tussen verschillende soorten 
deelljes die de reflecties veroorzaken en bij het uitvoeren van de maximum 
likelihood schatting de uitdrukking gecorrigeerd wordt met een product van 

15 respectieve tweede fistctoren voor de waarschjjnlijkheid dat niet tevens een 
reflectie met een dominerende waarde van de e:^nschap door 
respecfcievelijke soorten deeUjes, die de gemeten waarde zou maskeren, deel 
uitmaakt van de reflecfciemeting. 

7. Inrichting voor het meten van parameters die een 

20 deeltjesgrootteverdeUng van deelljes in een vloeistof beschrijveii, welke 
inrichtii^ omvat 

- een vloeistofkanaal; 

- middelen om een ultrasone bundel op te wekken in het vloeistofkanaal; 
* middelen om een eigenschap van een reflectie van de bundel aan een 

26 deeltje in het vloeistofkanaal te meten; 

- een verwerkingseenheid ingexidit voor het uitvoeren van een m aximum 
likelihood schatting van de parameters, gezien een combinatie van de 
gemeten waardes, op basis van een uitdrukking voor een waarschijnlijkheid 
van de gemeten waardes als fdnctie van de gemeten waardes, welke 

30 uitdrukking een eerste factor bevat voor de waarschijnlijkheid van een 
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reflectiemetixig waarvan een refleciie met de gemeten waarde onderdeel van 
uitniaakt, geoorri^^erd met een tweede £a.ctor voor de waarschqnlijkheid dat 
niet tevens een refleciie met een dominerende waarde van de e^;enschap, 
die de gemeten waarde zou maskeren, deel uitmaakt van de reflectiemeting. 
5 8. Computer prc^rammaprodncfc met instructies voor te^ 

parameters die een deellgeagrootteverdeling van deellges in een vloeistof 
beschrijven, aan de hand van van een reeks reflectiemetine^n, in elk 
waarvan een signaalbimdel in de vloeistof wordt opgewekt en een waarde 
(A) van een eigenschap van een reflectie aan een deeltje in de signaalbimdel 

10 wordt gemeten; welke instructies zijn ingericht voor het nitvoeren van een 
maximum likelihood schatting van de parameters gezien een combinatie 
van de gemeten waardes (A), op basis van een uitdrukking voor een 
waarschijnlijkheid van de gemeten waardes als functie van de gemeten 
waardes, welke uitdrukking een eerste £Eictor (P) bevat voor de 

15 waarsdignlijkheid van een refLectiemeting waarvan een reflectie met de 
gemeten waarde onderdeel van uitmaakt, gecorrigeerd met een tweede 
&ctor (Q) voor de waarschijnlijkheid dat niet tevens een reflectie met een 
dominerende waarde van de eigenschap, die de gemeten waarde zou 
maskeren, deel uitmaakt van de reflectiemeting. 

20 
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